New theory of the 7-a phase transition in Ce: quadrupolar ordering 
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We present a theoretical model of the "isostructural" 7 — a phase transition in Ce which is based on 
quadrupolar interactions due to coupled charge density fluctuations of 4/ electrons and of conduction 
electrons. Conduction electrons are treated in tight-binding approximation. The 7 — a transition is 
described as an orientational ordering of quadrupolar electronic densities in a Pa3 structure. The 
quadrupolar order of the conduction electron densities is complementary to the quadrupolar order 
of 4/ electron densities. The inclusion of conduction electrons leads to an increase of the lattice 
contraction at the 7 — a transition in comparison to the sole effect of 4/ electrons. We calculate 
the Bragg scattering law and suggest synchrotron radiation experiments in order to check the Pa3 
structure. The theory is capable of accounting for transitions to phases of non-cubic symmetry, but 
it is not sufficient to describe the magnetic phenomena which we ascribe to the Kondo mechanism. 
We also present a microscopic derivation of multipolar interactions and discuss the crystal field of 
7-Ce. 



I. INTRODUCTION 

The understanding of the nature of the a phase and 
the apparently isostructural transition between the cubic 
7- and a- phases in cerium is a long-standing problem 
Numerous experimental data present an outstand- 
ing challenge for the explanation by theoretical models 
[p|-|ri"|| . The most important question to be answered by 
theory is to find the driving force of the "isostructural" 
7 — a transition and to explain also the existence of the 
other non-cubic phases of Ce. 

The theory should also describe the change of magnetic 
properties of Ce at the 7 — a transition Jl^|l| which are 
reminiscent of an insulator-metal transition of 4/ elec- 
trons ||. Recently the Mott transition has been recon- 
sidered by several authors using electronic band structure 
calculations with thermodynamic extensions (see 
for a review Ref. Q). There, the a phase is described as 
a regular band state formed by 6s, 5c? and 4/ electronic 
states while in 7-Ce different degrees of localization of 4/ 
states are suggested and investigated. Two face centered 
cubic (fee) phases of cerium are attributed to two local 
minima of free energy which develop for the same crystal 
structure (Fm3m). 

Assuming the localized nature of the 4/ electrons 
throughout the 7 — a phase transition, one can under- 
stand the magnetic properties on the basis of two princi- 
ples. These are the singlet ground state of one 4/ electron 
and the energy gap (characterized by the Kondo temper- 
ature Tk) separating the ground state from a manifold 
of excited magnetic states. So far these properties have 
been treated JlO| , |TT]] on the basis of the Anderson impu- 
rity Hamiltonian |13| which implies an antiferromagnetic 
Kondo spin interaction. In order to describe the vol- 
ume contraction |pj at the 7 — a transition, the theoret- 



ical models ||[||7]] exploit the volume dependence of the 
Kondo temperature. Such a "volume collapse" leads to a 
phase instability without symmetry change and is inter- 
preted as an isostructural transition. However within the 
Kondo theory scenario, the existence of the other, non- 
cubic phases of Ce, remains unexplained. In the last few 
years the validity of this approach has been questioned by 
photoemission spectroscopy experiments JI^18[| where 
the predicted temperature dependence of the intensity 
of the Kondo resonance at the Fermi level has not been 
observed. Recently it has appeared that the Kondo vol- 
ume collapse model can not be applied to YbInCu4 p~9| ] 
which exhibits a 0.5% volume expansion during another 
isostructural phase transition (20| though Yb is the f-hole 
analogue of Ce. In addition, taking into account the ther- 
mal expansion of YbInCu4 above the phase transition 
temperature T = 42 K one concludes that the Kondo 
temperature is not a unique function of cell volume [fl9| . 

An alternative theoretical model of the 7 — a transi- 
tion has been recently proposed by the present authors 
|filf . The theory is capable of accounting for transitions 
to phases of non-cubic symmetry. The quantum mechan- 
ical electric quadrupole interaction between 4/ electrons 
on the fee lattice is proposed as the driving mechanism of 
a phase transition. The 7-phase is characterized by the 
absence of spatial orientational order of the quadrupo- 
lar densities, the space group is FmZm. In the a-phase 
the quadrupolar densities order in a Pa3 structure. No- 
tice that this change from Fm3m to Pa3 conserves the 
fee structure of the atomic center of mass points and is 
solely due to orientational order of the quadrupoles. This 
phase transition is accompanied by a contraction of the 
fee lattice, however the theoretical estimation of these ef- 
fects in ^lj (in the following we will denote this reference 
by I), is about an order of magnitude smaller than the 
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experimental result Q. In addition the treatment of I 
does not indicate the existence of a critical end point of 
the phase separation line (7 — a) in the P — T (pressure- 
temperature) phase diagram. 

In the present paper, we will extend the theoretical 
model of I by taking into account the polarization of 
(6s5d) 3 conduction band electrons. The polarization can 
be considered as a screening process of the quadrupolar 
density orientations of the 4/ electrons and results in turn 
in a complementary ordered Pa3 structure build up from 
conduction electron quadrupolar densities. The conduc- 
tion electrons will be described within the formalism of 
tight binding approximation. 

The paper comprises the following sections. We start 
(Sect. 2) with reconsidering 4/ electrons and extending 
the treatment of I by taking into account the radial de- 
pendence of the 4/ electron density. Next in section 3 
we describe the conduction electrons in second quantiza- 
tion with basis functions in tight-binding approximation. 
We derive the multipolar interactions among conduction 
electrons as well as interactions with 4/ electrons. In 
considering conduction electrons we have to distinguish 
between on- site and inter- site interactions. Section 4 is 
devoted to a study of the crystal field which acts on the 
individual 4/ electron. It is found that the refinement 
by the radial density dependence of 4/ electrons does 
not improve but rather spoils the agreement between ex- 
periment and theory. On the other hand the inclusion 
of conduction electrons improves in turn the situation. 
In Sect. 5 we study the quadrupolar ordering of the cou- 
pled system of 4/ electrons and conduction electrons. We 
conclude that a state of lower free energy (in comparison 
with the disordered state) can be achieved by a comple- 
mentary ordering of 4/ electron and conduction electron 
quadrupolar densities. On a same atomic lattice site, 
high charge density regions of the 4/ electron correspond 
to low charge density of the conduction electrons and vice 
versa. Finally we present numerical estimates that the 
inclusion of the conduction electron quadrupolar order 
improves the magnitude of the lattice contraction at the 
Fm3m — > Pa3> transition. In order to propose an un- 
ambiguous experimental proof of the present theoretical 
model, we calculate the Bragg scattering law and suggest 
synchrotron radiation experiments for the Pa3 structure 
(These calculations will be presented in the full original 
version to be published in Eur. Phys. J. B, 2000). In 
Sect. 6 (Discussion) we recall the salient features of the 
present theory and situate it with respect to the con- 
ventional approaches that are based on the Kondo the- 
ory concepts. While we do not adopt the Kondo volume 
collapse models, we discuss the relevance of the Friedel- 
Anderson hybridization mechanism for the explanation 
of the magnetic anomalies in Ce and suggest a link with 
our theory of the electronic charge degrees of freedom 
driven structural phase transition. 



II. RADIAL DEPENDENCE OF 4F ELECTRONS 

In our previous paper ]2l| ] the quadrupolar coupling 
between 4/ electrons has been calculated by assuming 
that the electron on each lattice site n is localized on 
a sphere with a fixed radius r/ = 1.378 a.u. Here we 
want to extend the previous calculation by taking into 
account the radial dependence of the 4/ electron wave 
functions. Such an extension is necessary if we want to 
study the interaction with the conduction electrons (see 
next sections). Furthermore it is useful in assessing the 
validity of our previous calculations. 

We consider a face centered cubic crystal of N Ce 
atoms. Each atomic core possesses one 4/ electron. In 
the 7-phase the 4/ electron densities are orientationally 
disordered. The space group of the crystal is Fmim. The 
Coulomb interaction between two 4/ electrons (charge 
|e| = l) at positions R(n) and R'(n') near the lattice sites 



and ft' reads 



V{R(n),R'(n')) = 



1 



\R(n) - R'{r 



The position vector R(n) is given by 
R(n) = X(n)+r(n). 



(2.1) 



(2.2) 



Here X(n) is the lattice vector which specifies the cen- 
ters of the atoms on a rigid fee lattice, while r(n) is the 
radius vector of the 4/ electron; in spherical coordinates 
r(n) = (r(H),Cl(fi)), where f2 = (Q,<p). We perform a 
multipole expansion of V by using site symmetry adapted 
functions (SAF's) |25| which transform as irreducible rep- 
resentations of the cubic site point group Oh- 



V(R(n),R'(n')) = £ v AA , (n,n'; r,r') S A (n) S A ,(n'), 

(2.3a) 



AA' 



where 
v AA ,(ri,ri'; r,r') 



dfl{n) / dQ(n') 



S A (n)S A >(n') 



\R{n) 



R'(n')\ 
(2.3b) 



The SAF's S A (h), h = fl{n) are linear combinations of 
spherical harmonics Y™ |Pq| . The index A stands for 
(/,t), with t = (r, /1, k). Here I accounts for the angu- 
lar dependence of the multipolar expansion, T denotes 
an irreducible representation (in the present case of the 
group Oh), fJ- labels the representations that occur more 
than once and k denotes the rows of a given represen- 
tation. Expansion (2.3b) still depends on the instanta- 
neous radii r(n) and r'(n'). In Ref. I we have written the 
Coulomb interaction in the space of orientational state 
vectors \i), i = 1 — 14, of the crystal field. The wave func- 
tions (n\i) were taken as linear combinations of spin or- 
bitals Y™(Q)u s (s z ) m = —3, +3, where u s is the spin 
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FIG. 1. Calculated radial density distribution for 4/, 5d 
and 6s electrons of an isolated Ce atom. Valence density is 
superposition of (6s) 2 and 5d; r c is the close contact radius 
of 7-Ce. 



function, with s = ± for the spin projections s z = ±1/2 
on the z-axis respectively. (The consideration of spin or- 
bitals is necessary since we calculate the eigenvalues Ei 
of the cubic crystal field in presence of spin-orbit cou- 
pling.) Since at present we take into account the radial 
dependence of the orbitals, we consider basis functions 



(n,r\i) = Kf(r(n))(n\i), 



(2.4) 



where we have assumed that the function lZf(r) is the 
same for all i. The real radial function TZf(r) is obtained 
from a DFT (density functional theory) calculation of an 
atom of Ce within LDA (local density approximation) for 
J = 5/2 states. In Fig. 1 we plot the radial density for 
the outer electrons. For a non-relativistic hydrogen-like 
atom IZf(r) would correspond to the Laguerre function 
TZ n=i i =3 (r). 

The matrix elements of the interaction ( |2.3a ) are ob- 
tained as 



^'\n'V{R{n) : R'(n'))\j) il \f)r, 

AA' 



(2.5) 



where 



v AA , (n - n ) 



dr r 

2 , 



dr' r" 

2 , 



y.K){r)K}{r')v KK ,{n,n'- r,r') (2.6) 
accounts for the average radial dependence and where 



c a(u) = dn(i\h)S A (n)(n\j). 



(2.7) 



We use the superscript F in order to indicate that we 
have transitions between two 4/ states, i.e. F = (/,/)■ 
By summing V(R(n), R' (n 1 )) over all pairs of lattice sites 
n, n' ', the total Coulomb interaction operator is then ob- 
tained as 



Uff 



where 



ftrt' AA' 



F / — \ FF t ^ 

Pa W %A' ( n 



n')p£(rf), (2- 



Pi 



(«) = E c a(^)I^™^I™- 



Introducing Fourier transforms 



n 

where q is the wave vector, we get 

q AA' 



(2.9) 

(2.10a) 
(2.10b) 

(2.11) 



The multipolar interaction (2.8) or equivalently (2.11) 
can be separated into two parts. We recall that A = 
(Z,r). 

Firstly we consider the case where I ^ and I' ^ 0. 
In I it has been shown that some of the coefficients 
c A (ij) = cj(ij) for I = 2 and r = (T 2g ,fi = 1, k = 1 — 3) 
are different from zero. (We recall that i and j refer to 
quantum states of the 4/ electron) . Therefrom we have 
inferred the existence of quadrupolar (I — 2) density fluc- 
tuations caused by transitions among 4/ electron states. 
In labeling the quadrupolar T 2g functions, we recall that 
the functions S\ = S(i = 2,T 2g ,k) are proportional to the 
Cartesian components yz, zx and xy for k = 1, 2 and 3 
respectively. In the basis of real spherical harmonics Y" z °, 
Y™ c and Y, m ' s of Ref. @ (see also Eqs. (2.1a)-(2.1c) 
of I), these functions correspond to Y 2 ' s , Y 2 ' c and Y 2 2 ' s . 



Writing only the index k for A = (I 
the quadrupolar density operator by 



pI 



2,k), we denote 



(2.12) 



The interaction between quadrupolar 4/ electron densi- 
ties becomes 

U QQ = ll2J2Pk(^^(Q)Pk'(q)- (2-13) 

q kk' 

The explicit form of v££(q) is discussed in Appendix A. 
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Secondly we have the case where I ^ and V = or 
vice versa. This means that we consider a multipole I 
on a given lattice site while the surrounding multipoles 
on neighboring fee lattice sites are taken in spherical ap- 
proximation. This interaction contributes to the crystal 
field. The crystal field has unit cubic symmetry, the low- 
est nonzero value of I is 4 and r = {A\ g , 1), where A\ g is 
the unit representation of the cubic site group Oh- The 
crystal field contribution from 4/ electrons at site n is 
then given by 



v£ F (f£) 



12 



V Iai 9 



F 

0,Ai s 



F 
Pl.Au 



(n) 



(2.14) 



Here we have restricted ourselves to the 12 nearest neigh- 
bors on the fee lattice, and I = 4,6, ... . 

The leading contributions to XJ*' are then given by 



u, 



ff 



CF' 



(2.15) 



where U f J F = V c 



Q F {n). In the following of this sec- 



tion we will discuss the physical consequences of the 



term UqL, expression (2.13), also called orientational 



QQ' 



pair quadrupolar interaction. We will give a discussion 
of the crystal field in Sect. 4. 

Previously (I) it has been found that the quadrupole- 
quadrupole interaction matrix v FF (q) becomes diagonal 
and has a twofold degenerate negative eigenvalue called 
X x + at the X point of Brillouin zone (BZ) of the fee 
lattice (see also Appendix A). This attractive interaction 
induces an orientational ordering of the quadrupolar den- 
sities in a Pai structure. A condensation scheme for the 
phase transition Fm3m — > Pa3 reads: 



Pi (C ) = Pi K ) = f>2 (C )=p t VN?0, (2.16a) 

pZffi) = pii€) = ft iff) = °- ( 2 - 16b ) 

Here the bar stands for a thermal expectation value, 
while p F is the order parameter amplitude. The above 
condensation scheme corresponds to one of eight possible 
domains of Pa3. The wave vectors q£ , q x and q x stand 
for (2tt/o) (1,0,0), (2tt/o) (0,1,0) and (27r/a)(0, 0, 1) re- 
spectively, where a is the cubic lattice constant. For 
each arm of the star of *q x — {q x ,q x ,q x } there are 
two basis functions p F . Hence the functions p F (q x ) of 
the condensation scheme (2.16a,b) form a basis of the six 



dimensional irreducible representation X§ of the space 
group Fm3m. In real space the ordering implies four 
sublattices of simple cubic structure as shown in figure 3 
of Ref. I. 

In I, where we have taken a fixed radius r/ for the 
4/ radial distribution, we found the eigenvalue X x + = 
—3491 K (Kelvin) and a phase transition temperature 
T\ = 85.6 K. At present we have calculated the eigen- 
value X x + = —4:'-f FF (see Appendix A) with the radial 




RADIUS (a.u.) 

FIG. 2. Radial dependence of quadrupole-quadrupole in- 
teractions uaa' (n, ft'; r, r'), 
where r = r' , and A = A' = (I = 2,T2 9 ,fc); ft = ft' (intra-, 
k — 1, 2, 3) and ft 7^ ft' (inter- on fee lattice) with ft = (0, 0, 0), 
ft' = (0/2X0,1,1), k = 1. 



DFT calculation with LDA. Although the atomic 4/ elec- 
tron density is small beyond the close contact radius 
r c = a/(2y/2) = 3.448 a.u., the opposite holds for the 
inter- site interaction potential VAA'(n,n'; r,r') and its 
first derivative (see below), which increase substantially 
when r — ► r c and r' — > r c (see Fig. 2). We have inves- 
tigated several models for the radial integral (2.6). In 
model 1 we consider the radial integration in the range 
< r < r c , that is without any overlap of the atomic 4/ 
electronic densities of neighboring atoms, and we obtain 
A 



-2121 K, next (model 2) we have allowed for an 
overlap between neighboring sites and extended the in- 
tegration over the range < r < 00. We then obtain 
= —2478 K. Finally (model 3) we assume again the 
integration range < r < r c but renormalize the 4/ elec- 
tronic density to unity. The result for X x + is -2682 K. 
Comparing these values with X x + = —3491 K obtained 
for the calculation with the characteristic radius rj, we 
conclude that a refinement of the theory in smearing out 
the radial extension of the 4/ electron distribution does 
not increase the strength of the quadrupole-quadrupole 
interaction and consequently does not increase the tran- 
sition temperature. 

So far we have considered multipolar interactions on a 
rigid lattice. In order to account for the lattice contrac- 
tion at the 7 — a phase transition, we have to include 
lattice displacements. In I we have shown that the inter- 
site quadrupolar interaction is modified by lattice dis- 
placements u v (n) . The correction to the potential reads 



dependence lZf(r) in Eq. (2.6) taken from an atomic 



Uqqt 



\Y1 J2J2 v lkk'(n-n';r,r') 
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x5 2 fc (n) Sl'(n<) [u u (n) - u u {ft')) . (2.17) v' v FF %p) = {Nm)~ 1 ^ 2 A FF q v a cos (^) cos (M) 



Here S% are SAF's with Z = 2, T 2s , Jfe = 1 - 3. The 
coupling coefficients v' v kk , are given by the derivative of 
the quadrupole-quadrupole interaction with respect to 
lattice displacements: 



d 1 



dX v {n) \R(ji) _ R'(j 



(2.18) 



Previously (I) this expression was considered for 4/ elec- 
trons on a shell with the (fixed) characteristic radius 
r = r' = rf. Defining 

Vl(R{n),R'(n')) = E v' v kk , (ft - ft'; r, r') S k 2 \n) S k 2 ' \n') 

kk' 

(2.19) 

we now consider matrix elements with basis functions 



(2.4) and obtain 



(i\H^%>Vl(R(n),R'(n'))\j) H \j%, 



kk' 



where 



v'M{ft -n') = j drr 2 J dr< r' 2 K 2 (r) TZ 2 (r') 

xvUk'in-ft'; r,r'). (2.21) 

The correction to the quadrupolar interaction between 
4/ electrons then becomes 

U QQT = lj2'pk(n)v' v ^(fi-n')p^(H') 

nn' 

x [u v {ft) - u u (ft')] , (2.22) 

with summation over repeated indices k, k', v. We intro- 
duce the Fourier expansion 



i v (ft) = (Nm)- 1 / 2 J2MQ)^ 



X(n) 



(2.23) 



where m is the Ce mas s. Using definition ( 2.10a ), we 
rewrite expression (2.22) in Fourier space. In the long 
wavelength limit q — > and taking p close to the star of 
p , we obtain 

U QQT = * E Y!<kk F %P) Pk (P) f Pk (?) MO)- (2-24) 

pq v(k) 

Here the sum J^' refers to i/ = x, y for fc=3, to v = z,x 
for k=2 and to ^ = z, ?/ for fc=l. The coupling matrix is 
obtained as 



(2.25) 

with A FF — v' u FF (n — ft'), where v — x or y and X(ft) — 
X(n') = (a/2)(l, 1, 0) on the fee lattice (for more details, 
see Ref. I). The other elements of v' u FF (q,p) follow by 
symmetry considerations and permutation of indices. 

We consider expression (2.24) in the Pai ordered phase 
by using the condensation scheme ( |2.16a ). The lattice 
displacements are taken in the long wavelength limit 
where they are related to the homogeneous strains. Sym- 
metry implies that only longitudinal strains occur: 



lim iq y u v (q) 

<f->0 

Then Uqqt becomes 



V mNe vl/ , v = x,y,z. 



ff _ 



jy QQT 



-2aA FF (p F ) 2 J2^, 



(2.26) 



(2.27) 



which corresponds to a coupling of ordered quadrupoles 
(quadratic) to longitudinal strains. Notice that the sign 
of A FF is negative, this is a consequence of the repul- 
sive nature of quadrupole-quadrupole interaction. The 
strains give rise to an elastic energy of the cubic lattice: 

1 a 3 
Utt= [c° u (e 2 xx + e 2 yy + e 2 zz ) 



N 



^^^Vl^-xx^yy ~t~ tyy^zz ^zz^xx)\ • (2.28) 



Here and c° 2 are the bare elastic constants. The in- 
terplay of quadrupolar order and lattice displacements 
follows from the interaction Hamiltonian 



U=U TT + Uqq + Uqqj,. 



(2.29) 



Minimizing U[p F , e vv \ with respect to the strains for a 
given Pa3 ordered structure, we obtain e w , 



txx = £ yy = e zz = -8a 2 |A FF | K L (// ) 



(2.30) 



while the change of the lattice constant is given by 
Aa = e xx a. Here kl = {c\\ + 2c < ( 2 )~ 1 is the linear 
compressibility. Hence the present theory leads unam- 
biguously to a lattice contraction. In order to provide 
a numerical estimate of the lattice contraction, one has 
to calculate A FF , to estimate Kl from experimental re- 
sults and to calculate by the methods of statistical me- 
chanics the discontinuity of the order parameter at the 
7 — a phase transition. In I (fixed radius rf) we did ob- 
tain A FF = -445 K/a.u. (-841 K/A) Q. Now, for the 
models 1), 2), and 3) of the spatial radial integrals we ob- 
tain A FF =-498, -581 and -659 K/A respectively. We 
conclude that the present refinement in calculating A FF 
is not helpful in view of obtaining a larger value of A FF 
and hence of the lattice contraction Aa. Since our theo- 
retical value of Aa is more than one order of magnitude 
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too small to account for the ~ 15% volume contraction 
of cerium, we conclude that a significant process has so 
far been omitted in our treatment. In the following sec- 
tion we will study the effect of conduction electrons on 
quadrupolar interactions. 



III. MULTIPOLAR INTERACTIONS 

Here we will investigate about the existence of multi- 
polar interactions between the localized 4/ electrons and 
the (6s5d) 3 conduction electrons U~f c as well as multipo- 
lar interactions among conduction electrons U cc . 

As a consequence of Bloch's theorem the conduction 
electron states are classified according to the irreducible 
representations of the translational symmetry group of 
the crystal. Nevertheless, in the proximity of the nu- 
clei the corresponding wave functions can be expanded 
in terms of spherical harmonics. This fact reflects the 
importance of Coulomb singularities associated with the 
nuclei. In the following we will focus on the interac- 
tions inside the "muffin-tin" or touching spheres centered 
on the nuclei and for the description of itinerant states 
adopt the tight-binding approximation. In absence of a 
static magnetic field, conduction electronic states with 
spin projections s z — ±1/2 are degenerate. Hence we 
will omit the spin dependence of the wave function. The 
wave function of a conduction electron with wave vector k 
and band index a is then written as a linear combination 
of local atomic wave functions 

(R\k,a)=^ a (R) 

= 7y= E emR ' } E 7»m(*, a) cj> lm {R - X(n% (3.1) 

n' Ira 



where the position vector R is given by Eq. ( |2.2| ). 
The atomic wave functions are given by <j>i m (f) = 
Hi(r)Y™(£l). The expansion coefficients 7/ m (fc, a) and 
corresponding eigenvalues E{k,a) are obtained by solv- 
ing the secular equation (see e.g. |Q) 

E V $ - E & a ) S * V ®) 7V (k, a) = 0. (3.2) 



Here H(k) and S(k) are matrices of single particle Hamil- 
tonian and overlap, respectively. Here and in the follow- 
ing A = (l,m), 6 xx' — Su'S mm i. The eigenvalues E(k, a) 
refer to the energy band spectrum of the conduction 
electrons. Notice that the wave function a (R) sat- 
isfies the Bloch condition. In the following we shall use 
an idealized basis set of orthogonal Wannier functions 
4>\(f{n)) = (R\\)fz without overlap, 



n(A|A'),t' — SxX'Sfin'- 

The coefficients 7a satisfy the relation 



(3.3) 



E 7A ( fc i a )7A (&',«') =5%%,5aa! 



The functions ipt (R) are then normalized: 



(k, a\k', a') = SgpSaa'. 



(3.4) 



(3.5) 



The corresponding electronic Hamiltonian in second 
quantization is 



ka 



(3.6) 



where at and ar are creation and annihilation Opera- 
te ka *" 

tors for one electron in state (k, a). 

In case of cerium, we will restrict ourselves to conduc- 
tion electrons of d type (I — 2, to=-2,...,+2, with prin- 
cipal quantum number n=5) and of s type (1=0, m=0, 
with n=6). Fig. 1 shows the radial density of the va- 
lence electrons of a Ce atom, it should also be a guide of 
the relative spatial extension of the localized 4/ and the 
conduction electrons in the crystal. 

Although the derivation of multipolar interactions pre- 
sented later in this Section is general in the following we 
imply the high temperature 7 phase of cerium. There the 
total electronic density associated with the conduction 
electrons has the symmetry of the space group Fm3m 
pl[ . At each lattice site the band electronic density 
has on average the unit symmetry (A\ g ) of the cubic 
point group Oh- It is important to realize that the on- 
site quadrupolar charge density fluctuations of T^ g and 
Eg symmetry do not interact with the totally symmetric 
charge distribution {A\ a ) rep resenting the ground state 
of cerium (see also Eq. j3.ll )). On one hand this allows 
us to use the Landau concept of Fermi liquid and to con- 
sider the conduction electrons described by Eq. (3.1) as 
quasiparticles where the ground state on-site interactions 
have already been taken into account through the coeffi- 
cients Ji m (k, a). On the other hand we can focus only on 
the quadrupolar interactions where conduction electrons 
are involved and consider these interactions separately 
from the ground state. 

Interactions involving conduction electrons will be de- 
scribed within the formalism of second quantization. We 
introduce field operators 



= J^a Ua {R\ka), 

ka 

&(R)=J2al a (ka\R). 



(3.7a) 
(3.7b) 



The operators W and ^' satisfy the usual anticommuta- 
tion relations for fermion fields. Although the derivation 
of multipolar expansions for these interactions is analo- 
gous, in particular as far as symmetry is concerned, to the 
procedure of Sect. 2, there is an essential difference. Since 
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the conduction electrons are not localized, they give rise 
to quadrupolar pair interactions on a same site (intra) 
and to interactions between different sites (inter). 

The Coulomb interaction between localized 4/ elec- 
trons at sites {ft} and the conduction electrons is given 
by 



dR' ^(R')(i\ n V(R(rt),R')\j).nt>(R'). (3. 



We observe that J dR' — > J2 n , I dr', where the integral 
extends over the volume of the cell ft! . In expression 
( |3.8[ ) we perform a multipole expansion of V similar to 
Eqs. ( 2.34 b) and then calculate matrix elements 



3 i(p-k)-X(n') 
N 

x ^E"a A' 2 (™ ~ " ) c a (ij) c A'hh ( k , P) P)- 

AA' hh 



(iy(k,a\V(R(n),R'(n'))\j) n \p,p) = 

F hh Z>\ „F 



(3.9) 



Here we define 

vl ft 2 (ft - n')= I drr 2 I dr'r" TZ 2 f (r) 



xuaa' (ft, ft'; r, r') TZ h (r) Ki 2 (r'), (3.10a) 



which is independent of the site, as is also the case for 
V AA'(^ ~ ™')- H ere r > = max(r,r'), r< = min(r,r') 
and <5aa' = Stt'^w ■ The inter-site coupling A i(ft — ft') 
still depends on the distance \X(ft) — X(ft')\, as follows 
from the translational invariance of the lattice. 

In additio n to t he multipole density of 4/ electrons 
p A (q), Eq. ( [2. 10a ), we define the multipole density of 
conduction electrons 



P L A 



{n) = ^Y,R L Me-^\ (3.12a) 

(q) = -^=^2^2at aaj :_^c AL (k,a;k~ q,P). 

(3.12b) 



a P k 



The interaction for the inter- site contribution of U^ c is 
then given by 



Jjfc 



Epa®^aa-(?>A'(9), (3T3a) 



with 



< L A> 



M0 = E e ™<i'(^)' 

The on-site part of U^ c is obtained as 



U fc \. t = C 

I intra 



F L 
A A 



(3.13b) 



(3.14a) 



mi 7712 

x CA'(Zimi, Z2TO2), (3.10b) 



where 



and 

c k ,{hm x , i 2 m 2 ) = J Y^*(n)s A/ (n)Y^(n)dn. 

(3.10c) 

In the following we will introduce a single index L for 
(hth), writing t>ff, for £, 2 , c A 'l for ca'Ws etc.. We 
notice that the index Zi (Z2) takes the values and 2 cor- 
responding to s and d electrons. We recall that c^(ij), 
referring to the 4/ electron transitions, is given by ex- 
pression (2.7). 

Since the conduction electrons are delocalized, we will 
have to distinguish between interactions where ft =/= ft' 
(inter-site), and where ft = ft' (on-site). In the first case 
vaa 1 (ft, ft'; r, r') is given by an expression of type (Otj), 
in the second case we have 

v AA ,(rt = ft'; r,r') = I ' dSl ( dO. 1 1 g A (0) S A >(Q.') 



An 



r V+ l ) 21 + 1 



Saa>, (3.11) 



Cl L A = vi L k (ft = ft'). 



(3.14b) 



Here we have also used the orthogonality relation ( [3.11 ) 
In expressions (3.13a) and (3.14a), summation is under- 



stood over indices A, A', L (l\, I2). 

In a similar way we treat the multipolar interactions 
between conduction electrons. We now start from the 
expression in operator representation 

U cc = i JdR' fdR&(R')&(R)V(R,R')V(R)y(R'). 

(3.15) 

We have to consider matrix elements 
(k,a\(k',a'\V(R,R')W,0)\P,P') = 
EEEE^a %(n- n') c AL (k, a; p, 13) 

nn' AA' L L' 

xc A , L ,(k',a'; p',P') 
where 



,i(p-k)X(n) A(p'-k')X(n') 



N 



N 



(3.16) 



Va a> in ~ ft') = I dr r 2 J dr' r' 2 K h (r)K h (r) 

x v AA i (ft, ft'; r, r') (r')T^ (r'). (3.17) 
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Again we distinguish inter- site and intra- site interac- 
tions. We obtain for the inter- site contribution 



U c 



1 \ t t 

> al a- ar, - 

2N ' ka k ' a ' k qt 



qf3 a k+q(3> 



kk' q 



x v a A' (?) °AL(k, a; k + q, 0) 
x c A > L>{k',a';k' - q,j3'), 



with 



V A A 



:{q) = Y,^ X ^vi%{h). 



(3.18a) 



(3.18b) 



In Eq. ( 3. 18a| ) summation is understood over the indices 
a..., Lfaflz)... , A... . Expression ( |3.18a ) can be rewrit- 
ten as 



(3.19) 



where pj^(q) is given by Eq. ( 3.12b| ) and where r) is the 
normally ordered product operator such that all the 
are placed to the left and all a to the right in the product. 
The on-site contribution is given by 



tjcc\ L 

U lintra ™ rj A A 



with 



Caa = "a a (" = "') 



L L' f-t 



(3.20a) 



(3.20b) 



So far the present formalism is general as far as multi- 
poles are concerned. In the following of this section we 
will study the interaction between quadrupoles (I = 2, 
I' = 2). In Sect. 4 we will study the crystal field (I = 4, 
i' = 0). 

We consider the three quadrupolar components of T2 g 
symmetry and write the index k for A = (/ = 2, k), 
with k = 1 — 3. The quadrupolar density of conduction 
electrons becomes 



a/3 p 



(3.21) 



We observe that Ckhip, ot;p — q, (3) comprises the factor 
Ck(hmi, ^2^2) which is defined (compare with expres- 
sions ( |3.10b| ) and fl3.10cj )) by 



c k (hmx, l 2 m 2 ) = J Y l ™ 1 %n)S£(n)Y£*(n)tKl. (3.22) 
The quadrupolar contributions from the interaction po- 



tentials (3.13a), (3.14a) between 4/ electrons and con- 
duction electrons then read 



TABLE I. Calculated coefficients Ck(l\_m ly £2^2), 
k = 1 — 3, li = h = 2; m,i (i — 1,2) stands for the in- 
dices (mi,c) or (m,i,s) of real spherical harmonics [25] Those 
functions which are not quoted here give zero contributions. 



(m) 


(m') 


r=(T 2g ,l) 


(T 29 ,2) 


(T 29 ,3) 


(0,c) 


(be) 





0.09011 





(0,c) 


(bs) 


0.09011 








(0,c) 


(2,s) 








-0.18022 


(l.c) 


(bs) 








0.15608 


(l,c) 


(2,c) 





0.15608 





(l,c) 


(2,s) 


0.15608 








(1,8) 


(2,c) 


-0.15608 








(bs) 


(2,s) 





0.15608 






u, 



fc 



= Y,Pk(tfvkk>(?)Pk>(q), (3.23a) 



and 



U, 



fc 



QQ 



= C**£/>*W /£($)■ (3.23b) 
The contribution of the quadrupolar pair interaction po- 



tential to expression ( |3.8| ) is then given by 



U, 



QQ 



QQ 



U, 



fc 

QQ 



(3.24) 



The quadrupolar interactions between conduction elec- 



trons follow from Eqs. (3.19) and (3.20a): 

JJCC I 



nter = \ll *® ^Pkitf P% '(5)) , (3.25a) 
9 

^0 lintra = \ C * ' * E vUitf Pk' ®) • t 3 ^) 



The contribution to expression (3.15) is given by 

JJCC _ JJCC I , JJCC I 

QQ - W QQ\ inter WW I intra 



(3.26) 



In order to select the contribution from conduction elec- 
trons, we have studied the coefficients Ck(hmi, limi). 
We find it convenient to use as basis functions for the 
conduction electron states the real spherical harmonics 
(see Sect. [nj). Then the s electron state is \l = 0, m = 0) 
while the five d electron states are |2, m = 0); |2, (to, s)), 
|2, (to, c)), m = 1, 2. There are no transitions between 6s 
states, the transitions between 6s and bd states are 



(0,0|S 2 1 |2 ) (l lS )) = (0,0|5 2 2 |2,(l, C )) 
= <0,0|S 2 3 |2,(2,s)) 



-in 



(3.27) 



and zero otherwise. The transition elements between 5d 
states are quoted in Table 1. 
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TABLE II. Interaction parameters C AB , X AB and A AB 
calculated with s— , d— , /— atomic radial distributions and 
lattice constant a=9.753 a.u.(7-Ce); A — (hh), B — (I1I2), 



F = 


(ff)- 








hh 


i'A 


C (in K) 


A (in K) 


A (in K/A) 


ff 


ff 




-2121 


-498 


sd 


ff 


±4408 


-3924 


-950 


dd 


ff 


75389 


-6025 


-1459 


dd 


dd 


94559 


-17695 


-4283 


ds 


ds 


21297 


-7547 


-1822 


ds 


dd 


±34858 


-11545 


-2792 



Since the coefficients Ck(hmi, I2ITI2) always occur in 
conjunction with interaction matrix elements v l k l2 l k } 1 2 
i v k k' ) ' we can immediately select the relevant matrix 
elements. We observe that the indices l\ and h in the 
interaction matrix elements then directly refer to s or d 



sd dd 
k k> I 



„// sd 
J k k' 



electrons and we will adopt the notation v ^ 
etc.. The indices /, s and d of these interaction matrix 
elements refer only to the radial dependence, they are 
irrelevant for symmetry considerations. The structure of 
the intersite quadrupole-quadrupole interaction matrices 
are investigated in Appendix A. From Eq. ( |3.11 ) it fol- 
lows that the on-site elements C k k , are diagonal in fc, 
k' and equal for k— 1,2,3; the same holds for C k k , . In 
the following we omit the indices fc, k' and write C^ sd , 
(jsdsd e ^ c Numerical values are given in Table 2. By 
introducing quadrupolar pair interactions on a same site 
we generalize the concept of spherically symmetric on- 
site electron-electron repulsions which is a characteristic 
property of strongly correlated electron systems. 

In order to treat the effect of conduction electrons on 
the lattice contraction, we have studied the coupling of 
quadrupole-quadrupole interactions with lattice displace- 
ments within the tight- bind ing approach. We start from 
expression Uqqt, Eq. ( 2.17 ), and remind that u v (n) and 
U u (n') refer to lattice displacements at different site s. W e 
first consider matrix elements of VJ, expression ( 2.19 ), 
between 4/ electron and con duct i on elec tron states. We 
proceed in analogy with Eqs. ( |3.8| )-( 3.13b ) but retain only 
inter- site contributions. We then find (compare with Eq. 
(5.13) of I) in the long wavelength limit q — > 0, and p close 
to the star of p x : 



Uqqt = ^Y,H k(q,P) P F ki-P) p L k (p) uM- 

pq v(k) 



(3.28) 



Here the sum refers to v = x, y for /c=3; to v = z, x 
for k = 2 and to v = z, y for k = 1. The coupling matrix 
is obtained as 



v » 33iq ' lf) = 7W^ A ^ acos {— ) coa {— ) 



(3.29) 



with A F L given by v'^^in — n'), where v = x or y, k = 3 
and X(n)-X(n') = (a/2)(l, 1,0) on the fee lattice. Here 
we have defined 

v' v l L k ,{n-n') = Jdrr 2 Jdr'r' 2 n){r) 

x v' v w (n - ff; r, r') K h (r')TZ h (/), (3.30) 



with v' ukk ,(n — n';r, r') given by Eq. ( |2.18| ). 

The matrix elements of V' between conduction elec- 



tron states are treated in analogy with Eqs. (3.15)-(3.19). 
Now only inter- site terms occur. In the limit q — > and 
p close to the star of p x we find 



U QQT = ^T,' v ^k(P,<l)v(p L k(-p)p L k\p)) uM- 

(3.31) 



pq iv(fc) 



Here again we have the same relatio n bet ween the indices 
v and k as was the case for Eq. ( [3.28j ). The coupling 
matrix reads 



1 



-.A 



fP x a\ (P v a\ 
q^acos j cos ) , 

(3.32) 



with A LL given by v' vk k (ft — n'), with v — x or y, for 
k = 3. We finally quote the definition 

<k k> (n "«') = / dr r 2 J dr' r' 2 K h (r)K h (r) 

x v' vh An-n';ry)nry)n v y), (3.33) 



where we use expression ( 2.18| ) . Numerical values of A L L 
are given in Table 2. 

Taking into account the contributions due to conduc- 
tion electrons, we see that the interaction Hamiltonian 
( |2.29 ) has to be replaced by 



U — Utt + Uqq + Uqqt, 



where 



Uqq = U f Q f Q 
Uqqt = U f QQT 



U, 



fc 

QQ 
U, 



u, 



fc 

QQT 



QQ' 
U, 



QQT ^ 



(3.34) 

(3.35a) 
(3.35b) 



Before studying the quadrupolar ordering and the ac- 
companying lattice contraction (Sect. 0) during the 7 — > 
a phase transition, we will next investigate the crystal 
field of 7-Ce in the presence of conduction electrons. 



IV. CRYSTAL FIELD OF 7-CE 

In the disordered 7 phase there are only charge den- 
sity fluctuations of quadrupolar type and the quadrupo- 
lar Hamiltonian (3.34) (or (2.2E)) averages to zero. The 
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first nontrivial orientational interaction then corresponds 
to a crystal field Hamiltonian. In Ref. I we have de- 
fined the crystal field of 7-Ce as the potential experi- 
enced by a single 4/ electron at a site ft when spher- 
ically symmetric contributions (l'=0) from nuclei, core 
electronic shells, conduction electrons, 4/ electrons at 
the twelve neighboring sites ft' on the fee lattice and sim- 
ilar terms from the homogeneous electronic density in 
the interstitial regions are taken into account. In I we 
have shown that crystal field effects are reduced to a sin- 
gle particle term to which we have added a spin-orbit 
coupling for the 4/ electron. We first generalize the re- 
sults of I by taking into account the radial de pen dence 
of the 4/ electron density. We start from Eq. ( [2.8|) with 
A = (I = 4, Ai g ) = Ai and A' = (I = 0, A lg ) = 0. Since 
we are dealing with a Coulomb potential and a spheri- 
cally symmetric charge distribution, the coupling func- 
tion v\ 1 o(n,n';r,r'), Eq. ( 2T3b| ), does not depend on r' 
and we observe that Eq. (2J3) can be written as 



where 



F Ft- -i\ F • r\ 

«Aj i n ~ n ) = ^AiO ' Qf, 



drr 2 TZ 2 f (r)vA l0 (n,n';r,r') 



(4.1a) 



(4.1b) 



is the same for all 12 neighbors. We obtain q < 0. 
The charge in units e of the 4/ electron at a neighboring 
site is given by 



Q f 



dr' r' 



'Tl}(r'). 



(4.1c) 



In our model Qf = 1. However, if we distinguish two 
regions in the crystal, the first inside muffin-tin (MT) 
spheres and the second in the interstices, then in the 
MT-region Qf < 1 . The other contributions to the crys- 
tal field are dealt with similarly. Since the interaction 
parameter vj[ ,0 remains the same it is only the charges 
Qi (i stands for core, nucleus, conduction electrons and 
interstitial contributions) which we shall take care of. 
Collecting the contributions from the various charges Qi 
together with Qf, we obtain for the crystal field at site 



where Ai = (I = i,Ai g ), 



Bf 



12 



r\ Ft 

QeffCV Al0 



and 



(4.2a) 



(4.2b) 



(4.2c) 



Here e refers to the electron charge at site n (e = —1) 
and Q e f f to the effective charge of the surrounding neigh- 
borhood. From electrostatic considerations we find (see 



Appendix A of I) that Q e ff = (1 + %int)QMT, where 
Q mt is the total charge inside a MT-sphere (which is al- 
ways positive) and where the factor x int ss 2.853 accounts 
for the charge contributions from interstices for touching 



MT-spheres. The coefficients cj^ , Eq. (4.2c), are quoted 



in Appendix A of I, they are diagonal in the basis (A.9- 
A.ll) of I. The expressions (4.2a)-(4.2c) represent a re- 
finement of the crystal field calculations, Eqs. (A.3)-(A.6) 
of I, where now a radial distribution of the 4/ electron is 
taken into account. 

In I we have calculated the crystal field coefficient B? 
(we used the notation A for B* ) for a fixed radius r / . The 
obtained value £?^=346 K corresponds to rf — 1.156 a.u. 
and not to Tf — 1.378 a.u. as was quoted erroneously 
in I. With a radial distribution as specified in models 
1), 2), and 3) of Sect. |rj, we obtain B? =970, 1403 and 
1104 K, respectively. This implies that the strength of 
the crystalline electric field of the 4/ electron should be 
2.8-4 times larger than calculated in Ref. I and there- 
fore by the same factor (2.8-4) larger than the crystal 
field measured experimentally |2^ , |24]| . However, expres- 
sion (4.2a) is implicitly based on the approximation that 
the distribution of conduction electrons at site n is spher- 
ically symmetric and hence the on-site interaction from 
/ = 4 multipoles is ignored. In the following we will show 
that the effective crystal field decreases if we take into 
account the average on-site cubic distribution (I = 4) of 
conduction electrons. 



From Eqs. (3.14a,b) for the on-site interactions we ob- 
serve that the 4/ electronic density p 1 [ i couples with 
the conduction electron density of the same symmetry 
(r = Ai g , 1 — 4). In real space 
A = A' = (I = 4, Ai g ) ee Ai reads 



Eq. (3.14a) with 



U fc 



CfciXXW* Pa». (4.3) 



We ob tain the crystal field potential due to p\ {n), Eq. 
( |3.12a ), by replacing this quantity by its thermal average: 



(PaI 2 ) = X n l c Ald d(k, a; k, a) 5 h2 S h2 . 



(4.4) 



is the Fermi distribution 



Here n2 <J fi0 8 a/3 = {aia^^ 
and h = 2, l 2 = 2 takes into account the fact that only 
d electrons contribute to the matrix element c\ 1 i t i 2 . We 
then write (p A ) = (p\ ) ■ The corresponding crystal field 
at site n reads 



where 



Vl F {n)=Bf c p F Ki {n), 



(4.5a) 



(4.5b) 



with C£ Ai > 0- Collecting the contributions (4.2a) and 



(4.5a) we obtain the effective single particle potential act- 
ing on the 4/ electron at site n: 
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(4.6) 



We recall that B^ > 0. The sign of B^ c depends on 
(p\ ). We next investigate other interactions that affect 

(pi)- 

We consider the crystal field acting on p^in). By 
studying the inter- site interactions we fi nd that the ma- 
trix elements uaa' (n, n'; r, r'), Eq. (2.3b ), are negligibly 
small for I = I' = 4. On the other hand terms with I = 4, 
I' = are significant and we retain the interactions which 
involve 



dd m 

v At 



dr r 2 TZ 2 d (r) v Al (n, n';r,r'), 



(4.7) 



where vj~ ' < (compare with Eq. ( 4.1b )). Here again 
v Ai o i s independent of r'. Proceeding in analogy with the 
calculation of Vq f , we calculate the field due to the sur- 
rounding 4/ electrons, core electrons, nuclei, conduction 
electrons and interstitial contributions at site ft, thereby 
obtaining 



VS F (n) = B c p c Ax (n), 
where (n) — pj™ (n) and 
12 



B c 



(4.8a) 



(4.8b) 



with B c = 6093 K > 0. The positive sign of the field 
B c implies that (p^), calculated with the crystal field 
Vp F (n), is negative. Hence B^ c , Eq. (4.5b), is negative, 



which leads to a reduction of the effective crystal field of 



a 4/ electron in Eq. (4.6). 

The previous considerations indicate that the inclusion 
of conduction electrons leads to a reduction of the crystal 
field experienced by a 4/ electron. However, a more rig- 
orous approach should start with the crystal field Hamil- 
tonian of 7-Ce, 



U = % + U CF + U ao . 



(4.9) 



Here U$ is the "bare" electronic term ( |3.6| ), which in- 
cludes the kinetic energy and the spherically symmetric 
part of the electronic potential of conduction electrons, 
Uso = Sn^so(^) stands for the spin-orbit couplings of 
localized 4/ electrons (see Appendix A of I for details). 
Ucf is the crystal field comprising 4/ electrons and con- 
duction electrons: 

U C F = J2 V cp(n), (4.10a) 

ft 

V CF (n) = Vi F {n) + V^ c F (n) + V£ F {n) + V^ F {n). 



(4.10b) 



Here 



V§Un) = Cifi d 1 <(n)Vlf(n) 



is the on-site I = 4, I' = 4 interaction between conduction 
electrons. In a mean-field approximation Ucf leads to 
self-consistent crystal field potentials Vq F (n) and Vq F (n) 
for conduction and 4/ electrons, respectively. A quanti- 
tative calculation of these effects is beyond the scope of 
the present work. 



V. QUADRUPOLAR ORDERING 



In Section [II we have found that the system of local- 



ized 4/ electrons and of conduction electrons are cou- 
pled by means of the intra- site an d inter- site quadrupo- 
lar potentials ( 3.231 ) and ( 3.23a ). Such a bilinear cou- 
pling suggests that an ordering of 4/ electron density 
quadrupoles should imply an ordering of conduction elec- 
tron quadrupoles and vice versa. In Appendix A we 
have investigated the wave vector dependence of the 
quadrupole-quadrupole interaction matrices and the con- 
dition for quadrupolar order. The present section is di- 
vided into two parts. In part A we study the interplay 
of quadrupolar ordering of 4/ electrons and conduction 
electrons on a rigid cubic lattice; in part B, where we 
consider a deformable lattice, we show that quadrupolar 
ordering implies a lattice contraction. 
A. We are investigating the possibility of a condensa- 
tion of quadrupolar densities in a Pa3 structure. At 
the A-point of the BZ, the matrices of inter- site in- 
teractions {v£k'( < l)) an d ( v k k' become diagonal in 
k, k' and have two degenerate negative eigenvalues. We 
then are led in analogy with the condensation scheme 
( |2.16a ,b) for 4/ electrons to suggest the condensation 
scheme Fmim — > Pa3 for the quadrupole densities of 
conduction electrons: 



p*{€ ) = pii<% ) = p2(€ ) = p 
p%(q*) = psiqy) = pi(q?) = o, 



L ^ / N^Q; (5.1a) 
(5.1b) 



where L — (sd), (ds), (dd). 

We disentangle the various contributions to Uqq, Eq. 
Q3.35a ) , in the ordered a phase where we assume a simul- 
taneous condensation of quadrupolar densities of 4/ elec- 
trons and conduction electr ons. Taking into accou nt the 
condensation scheme ( 2.16a ,b), we obtain from Eq. ( 2.15 ) 



1 



If 



(5.2) 



N QQ 2 

where \ FF stands for X x + , the twofold degenerate nega- 
tive eigenvalue of v FF (q x ) (see Appendix A). Similarly, 
using in addition the condensation scheme ( |5.1a ,b), we 
obtain from Eq. ( 3.23b| ) 



N 



U, 



fc 



3 [2C 



Fds p F p- 



ds 



c 



FddpF -ddl 



(5.3a) 



(4.11) and from Eq. ( |3~23a] ) 
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fc 



TT. 



3 [2A 



F ds -F -ds . % 

p p + A 



F dd -F -ddl 
P P \ 



(5.3b) 



where A F ds and A F dd are the twofold degenerate negative 
eigenvalues of the matrices v (q x ) and v Fdd (q x ), re- 



spectively (see Appendix A). Finally expressions fl3~25bp 
and (3.25a) lead to 



— U cc I 

j\j QQ I intra 



and 

— U cc I 

M QQ I inter 



- \AC sdad {p sd ) 2 + AC sd 



dd —sd —dd 



[4\ sdsd (p sd ) 2 



+ C dddd (p dd ) 2 \ , (5.4a) 



f~ 4l\ S< ^ dd -sd -dd 

+X dddd, p ddy^ ^ 4b ) 



We observe that all coefficients C LL are positive while 
X LL are negati ve (s ee Table 2). Obviously the inter- 
site interaction ( 5.4b[) fav ors quadrupolar order while the 
intra- site coupling ( 5.4a ) disfavors quadrupolar order of 
conduction electrons. 

The leading quadrupolar interaction in the ordered a 
phase, 



Uqq = U, 



ff 
QQ 



U, 



fc 
QQ 



U 



fc 
QQ 



U, 



QQ\ inte 



(5.5) 



QQ I intra 

is a quadratic form in p F , p sd and p dd , in particular 
the terms (5.3a) and (5.3b) represent a bilinear cou- 
pling between the localized 4/ quadrupolar density p F 
and the conduction electron quadrupolar densities p dd 
and p sd . Since in tight-binding the (6s5d) 3 electrons are 
hybridized, we introduce the total conduction electron 
density 



-dd 



-sd 



p- = fT- + 2p s 
Defining average interaction coefficients 
1 



C fc 



f^jF ds _|_ ^jf d<t\ 



Fdd\ 



Xf- = l(X Fds + X Fdd ), 



1 



2(C dsd8 + C 
o 



dd dd 



(jds dd^ 
\ ds dd\ 



6 

^cc ^_ ^ s _|_ X dd dd \ds dd^ 

3 

we approximate Uqq by an effective interaction 



(5.6) 

(5.7a) 
(5.7b) 
(5.7c) 
(5.7d) 



^U QQ » I [\ FF (p F f + 2A^p F f + A cc (p-) 2 ] , (5. 



where 



A fc = C fc + x fc 



(5.9a) 




FIG. 3. Pa3 structure of the ordered a phase with 4 
sublattices (nl — n4). Grey quadrupoles correspond to inner 
4/ electron densities, white quadrupoles - to outer conduction 
electron densities with the opposite sign as shown in figure 4. 




FIG. 4. Quadrupolar density distribution of the 4/ elec- 
tron (inside sphere of r=1.3 a.u.) and of one conduction elec- 
tron (spherical radius r=3.4 a.u.), on scale. 



and 



A cc = (jee + X c 



(5.9b) 



We observe that A^ c > 0, A cc > 0. For a fixed value of 



we minimize Uqq with respect to p c and obtain 



A fc 



(5.10) 



We see that quadrupolar order p F of 4/ electrons pro- 
duces (as a type of mirror image) a quadrupolar order 
of conduction electrons p c of opposite sign. Herewith we 
associate the pictorial representation of Fig. 3 and 4. 
Regions of an excess density (+) of the 4/ electron dis- 
tribution overlap with a depletion (— ) in density of the 
conduction electrons and vice versa (Fig. 4). Substitut- 
ing the right hand side of Eq. ( 5.10| ) into (5.8) we get 
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N 



Uqq 



\xn + {Afc)2 



A< 



(p 



F\2 



(5.11) 



where we use the fact that X FF < 0. The bilinear cou- 
pling between quadrupolar densities of 4/ electrons and 
conduction electrons leads to an increase of the attractive 
interaction between 4/ quadrupolar densities. 

The results of the last two sections can be summarized 
by suggesting an effective Hamiltonian in the disordered 
7 phase, 

H = U QQ + U , (5.12) 
where Uo is the crystal field Hamiltonian, fl4.S|) , and 

U QQ = U f Q f Q + U f Q c Q + U% Q . (5.13) 

Here 



+c.p. 



(5.14a) 



4^ /c EE^(^®. ( 5 - 14b ) 

q k 

U QQ = ^EE^®^®' ( 5 - 14c ) 

q k 

where pj(q ), k = 1 — 3, is the quadrupolar density oper- 
ator (2.10a) of 4/ electrons while p k (q) is the quadrupo- 
lar density of conduction electro ns (3.21 ) with pi — 
pf d + 2pl d . The first term, Eq. ( 5.14a| ), represents an 
effective attraction at the X point of the Brillouin zone. 
The other terms in ( 5.14a ) (c.p.) follow by cyclic per- 
mutation of indices, i.e. when q£ — » , k changes 
2 — > 3 and 3 — * 1, etc. The direct lattice is face center ed 
cubic. In that case the condensation schemes ( 2.16a ,b) 



and (5.1a,b) correspond to a structural phase transition 
where the space group changes from Fmim to Pa3. 

Starting from the interaction ( |5.8|) we construct a Lan- 
dau free energy in the condensed phase. We take p\ = p F 
and p2 = p c as components of a two dimensional vector 
p = (pi,p~2) and define the matrix 



J = 



-|A FF | Af c 

A_Sc A cc 



(5.15) 



Inspired from the theory of orientational order in molec- 
ular crystals we write 



Here J-q is the free energy in the disordered phase: 

To = -T\nTr[e- Uo/T ], (5.17) 



where U$ is the crystal field, Eq. (|4.9|). The quantity 



X 



(5.18) 



is the inverse susceptibility matrix, T is the temperature 
and 



(pp)o = 



<(/^)) 2 >o 









((p c k {n)?)o 



(5.19) 



Here (...)o are single particle thermal expectation values 
that have to be calculated by means of Uq, Eq. (|4.9|). Cu- 
bic symmetry implies (see Appendix B for details) that 
the three quadrupolar components are equal and for ex- 
pectation values of conduction electrons we obtain: 



1 



N 2 



c kL (h,a;p,/3)\ 5 LL >5 kk , 



(5.20) 



Here Sl l' = &i x i' ± Si 2 y . The contributions and J 7 ^ 
stand for the third and fourth order terms in p. Sym- 
metry of the or der pa ramete r com ponents (see conden- 
sation schemes ( 2.16a| ,b) and ( 5.1a ,b)) implies that there 
exists a non-zero third order cubic invariant J-^ and 
hence the transition Fmim — + Pa3 is of first order. As 
a first approximation we investigate the possibility of a 
second order phase transition. Neglecting J 7 ^ and , 
we minimize T with respect to p F and to p c . We ob- 
tain two coupled homogeneous equations from which it 
follows again that p F and p c have to be of opposite sign. 
The compatibility condition leads to the transition tem- 
perature 



Tr 



\ (jff ~ Jco + ^(Jff + Jcc) 2 + 4(J/ C ) 2 ) , (5.21) 



where 

Jff = 1A FF 1 ((pjT(n)) 2 )o, Joe = A c c {(p c k (n)) 2 )o, (5.22a) 
J fc = Af°y/ ! {{ P F k {ri)Y) {{pl{n)Y) a . (5.22b) 

We observe that Tc > Jff and the inclusion of con- 
duction electrons leads to an increase of transition tem- 
perature. The calculation of the first order transition 
temperature T\, (T\ > 7c), and of the accompanying 
discontinuities p c (Ti) and p F (Ti) in the coupled order 
parameters would require the evaluation of the higher 
order terms J 7 ' 3 ' and J 7 ^. Such an endeavor poses not 
only analytical but also numerical problems that are be- 
yond the scope of the present work. We therefore will 
adopt below (Sect. 6) an alternative point of view and 
suggest experiments in order to check qualitatively the 
predictions of the theory. 

B We next investigate the lattice strains at the tran- 
sition to the Pa3 phase. We recall the result Uqq T , 
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Eq. ( |2.27| ) , for the coupling of strains to 4/ electrons 
with quadrupolar order. The same procedure, as out- 
lined in Sect. 2, is now applied to derive the coupling of 
longitudinal strains to conduction electrons. Using the 
condensation schemes (2.16a,b) and (5.1a,b) and taking 
into account the definitions of interaction parameters A 
(Table 2), we obtain from Eqs. ( |3.2S| ) and ( |3 1.2SQ 

±U£„ = -4a [2A"W + A f W] £ e - 



Similarly we get from from Eqs. ( |3.31 ) and ( |3.32[ ) 



N U QQT 



-2a 



^p^ds ds ^-ds~^ _|_ ^p^ds dd -ds —dd 

+A dddd (p dd ) 2 ]Y / ^- (5-24) 



The coefficients A are quoted in Table 2. We define the 
total interaction potential of ordered quadrupolar elec- 
tron densities coupled to longitudinal strains by 



Uqqt - U f Q f QT + Uqq T + Uq Qt . 



(5.25) 



Here the rig ht ha nd side terms are given by Eqs. ( [2.27 ) 
( |5.23 ) and ( 5.24 ), respectively. In analogy with expres- 
sions ( 5.7b , d) we define 



A /c 



— (A F ds 
2 V 



_/j^dsds _j_ p^dddd 

3 



A 



ds dd\ 



(5.26a) 
(5.26b) 



and use again definition (5.6). Then we approximate 
Uqqt /N by an effective interaction 

±U QQT » ~2a[A FF {p F f + 2A^p F f + A cc (p c ) 2 ] 

X^e w , (5.27) 



where a is the cubic lattice constant and where kl is 
the bare linear compressibility. From the numerical val- 
ues of Table 2 we obtain A FF = -498, K fc = -1205, 
A cc = —2966 (units K/A). Since p F and p c are of op- 
posite sign, the fc contribution on the right hand side 



of expression (5.29) leads to an expansion of the lat- 
tice while the // and cc contributions lead to a con- 
traction. From Table 2 we also find tha t A fc = 30516 K, 
A cc = 24033 K, and hence from Eq. fl5.10| ) we deduce 



that 



IP1 > P 1 



Hence we conclude that inclusion of 



(5 23) t ne conduction electrons in expression (5.29) leads to an 



increase of the lattice contraction by a factor 4.4 in com- 
parison with the effect of the 4/ electrons. A numerical 



calculation of the lattice contraction As 



ae. 



at the 



first order phase transition would require the knowledge 
of the discontinuities p c (Ti) and p F (Ti). Experimentally 
jp one finds in the pressure (P) - temperature phase 
diagram of the 7 — a transition in Ce indications of a 
critical end point where the lattice contraction vanishes. 
Although the transition Fm3m — > Pet3 is always of first 
order and hence leads to finite discontinuities of p c (Ti) 
and p F (Ti), there exists the possibility that the relative 
importance of the coefficients A^ c and A cc , which have 
opposite sign, changes as a function of (P, T) and hence 
the expression within square brackets on the right hand 
side of Eq. ( |5.29 ) could vanish or even change sign with- 
out the requirement that p c (Ti) = and p F {T\) = 0. 
This opens up a possibility of a lattice expansion during 
the discussed first order phase transition. Recent experi- 
ments on YbInCu4 which exhibits an isostructural phase 
transition similar to the 7 — a change observed in cerium 
(see p5| for references) indicate a 0.5% volume expan- 
sion at the transition. Within the present theory such 
behavior can be understood if for Yb in YbInCu4 the fc 
term responsible f or exp ansion prevails over // and cc 
contributions, Eq. ( 5.2S ). 



VI. DISCUSSION 



This relation allows us to express the longitudinal strains 
as function of the order parameters. The contributions 
Utt and Uqqt lead to a supplementary term in the free 
energy (compare with Eq. ( 5.1 6| ) ) : 



1 1 1 

JjFQQT = JjUTT + JjUqqt. 



(5.28) 



(There are still additional contributions from thermal lat- 
tice vibrations, but these are irrelevant here.) For a given 
quadrupolar order, i.e. p F and p c fixed, we minimize 
Fqqt\p f ' , p c , € vv ] with respect to the strains e xx , e yy and 
e zz and obtain 



e 2 z = 8a 2 k l 

[A FF (p F ) 2 + 2A fc p F p c + A cc {p c ) 2 ] , (5.29) 



The present paper is an extension of our previous 
model of the 7 — a phase transition in Ce (Ref. [Bl| or 
I) . We predict that the 7 — a transition is accompanied 
by a symmetry change Fm3m — > Pa3 in the electronic 
structure. The idea of such a proposal is borrowed from 
the theory of molecular crystals, where Pa3 structures 
due to orientational order of molecular mass distribution 
are not unusual. For example, such crystal symmetry 
occurs in Na0 2 @, N 2 |(]] and in solid C 60 fl). The 
conventional characterization of the 7 — a transition in 
Ce as a phase transition "without change of symmetry" 
is based on several X-ray diffraction experiments It 
is possible that domain formation in the a-phase has pre- 
cluded an identification of this phase as a Pa3 structure. 
If new experiments are done, particular attention should 
be given to the possible coexistence of domains J32||, 
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In comparison with our previous work, Ref. I, we have 
extended the model in two respects. First, we have in- 
cluded quadrupolar interactions between 4/ and conduc- 
tion electrons and secondly, we have calculated the rel- 
evant parameters of interactions (Table 2 and Sect. 2) 
by using the radial dependences of valence electrons ob- 
tained from a DFT-LDA calculation of a cerium atom. 
To our knowledge, in the literature there exists no mi- 
croscopic derivation of multipolar interactions between 
conduction and localized electrons in solids although the 
concept of quadrupolar moment of a 4/ electron shell is 
well established p3 34|. Therefore in Sect. 3 we have pre- 
sented a detailed calculation of multipolar interactions 
treating band conduction electrons in second quantiza- 
tion with wave functions in tight-binding approximation. 
Since the intersite quadrupole-quadrupole interaction is 
short ranged and anisotropic, special attention has been 
given to lattice site symmetry. While we have restricted 
ourselves here (see in particular Appendix A) to the case 
of an fee lattice (7 — > a Ce) , our procedure is general and 
can easily be extended to other structures. 

Our results can be understood on the basis of the fol- 
lowing generalized Hamiltonian which we ascribe to the 
7 phase: 



:C 



(6.3b) 



H 7 — Uq + Uqq + Uqqt + Utt- 



(6.1) 



Here Uq is the crystal field Hamiltonian (4.9). The term 
Uqq represents the quadrupole-quadrupole interaction 
comprising c ontrib utions from 4/ and conduction elec- 
trons (Eq. ( 3.35a )). The quadrupolar interaction be- 
tween localized electrons is due only to intersite contri- 
butions on the fee lattice, see Eq. ( 2.13 ) in Fourier space. 
The presence of conduction electrons leads t o two types 
of contributions: inter-site terms, given by (3.23a) and 
( 3.25a ) fo r Uqq ancl Uqq i respectively, and on-site terms, 
given by ( 3.23b| ) and ( 3.25b| ), correspondingly. Collecting 



inter- site and on- site terms separately, we write 

U QQ = U QQ\i„ter + U QQ\intra ' ( 6 ' 2 ) 

where 

q k,k' ^ 



1 



+< U y (q) Pi WPw («) + 2 v * V (£> V (p L k (tfp% % 

(6.3a) 

and the summation over L(L') — (sd), (ds), (dd) is im- 
plied. The structure of the interaction matrices V (q), 
A{B) = F or L, is discussed in Appendix A, the relevant 
eigenvalues X AB at q x are quoted in Table 2. The intra- 
site contributions are given by 



where the calcula ted parameters C AB ar e quote d in Table 
2. Next in Eq. (6T), Uqqt, given by ( 3.35b ), is a cor- 
rection to Uqq for a deformable lattice, while Utt is the 
elastic energy of the cubic crystal in harmonic approxi- 
mation. In terms of homogeneous strains, Utt is given 
by Eq. ( |2.28| ) and Uqqt by Eq. (|J|). In principle a bi- 
linear coupling term Uqt between quadrupolar electronic 
and displacive degrees of freedom |35|] should be included 
in expression (7.1). This term which is known from Jahn- 
Teller phase transitions |36| can be essential if we want to 
describe transitions from a cubic phase with quadrupo- 
lar disorder to non cubic phases with ferro- quadrupo- 
lar order |37j. However for the transition to the Pa3> 
phase, which wc identify with a-Ce, the atomic center 
of mass positions still occupy a face centered cubic lat- 
tice and the term Uqt is irrelevant. Indeed, the driving 
force for the transition Fm3m — > Pa3 is the quadrupole- 
quadrupole interaction which becomes attractive (in re- 
ciprocal space) at the X point of the Brillouin zone. This 
fact leads to an orientational order of quadrupolar elec- 
tronic densities on four different sublattices (Figs. 3, 4). 
The term Uqqt then prompts a lattice contraction at the 
first order phase transition. Notice that the term Uqt is 
found to vanish for a wave vector q at the X point of the 
BZ. We insist on these facts since within our view, the 
electronic charge degrees of freedom, together with the 
lattice displacements, are the driving forces of structural 
phase transitions in Ce and related compounds. 

The Hamiltonian (7.1) is not sufficient to describe the 
magnetic phenomena that occur at the 7 — a transi- 
tion Jl2||l|. In accordance with Kramers' theorem the 
quadrupolar ordering, as described by Eq. (7.1) is not 
accompanied by a magnetic ordering. We then conclude 
that the addition of an Anderson Hamiltonian term H c f 
which takes into account the Friedel- Anderson hybridiza- 
tion between conduction electrons and 4/ electrons as 
well as the repulsive energy among the 4/ electrons on 
a same site is necessary (for a review see Q). Such 
a Hamiltonian leads to the disappearance of local mag- 
netic moments below a characteristic temperature Tk- 
The Kondo temperature Tk increases with increasing hy- 
bridization matrix element V . We then conclude that 
the lattice contraction accompanying the Pa3 quadrupo- 
lar ordering or the quadrupolar order enhances the hy- 
bridization V and hence increases Tk- This can lead to a 
situation where the structural 7 — * a transition and the 
demagnetization of the 4/ state occur at a same temper- 
ature T\, Notice however that within this scenario the 
process is driven by the structural (quadrupolar) transi- 
tion at T\ and not by the Kondo volume collapse |||^. 
In case where the enhancement of Tk is insufficient, the 
quadrupolar order and the concomitant lattice contrac- 
tion would occur at T\ without Kondo anomaly (disap- 



15 



pearance of magnetic moment). The condition Tk < T\ 
does not ensure that the Kondo anomaly actually occurs 
at lower T. Quadrupolar ordering has been observed 
in a number of Ce, Pr, Tm and U based compounds 
p3[ . A remarkable example is the magnetic semiconduc- 
tor TmTe with 4/ 13 ( 2 F 7 / 2 ) electronic configuration and 
Neel temperature 1^=0.43 K ||]. Although TmTe had 
been extensively studied before 1995, the phenomenon 
of quadrupole ordering below Tq=1.8 K was completely 
overlooked [^,^T). Numerous data on such compounds 
(CeAg 0, CeB 6 @, TmTe @, TmAu 2 @, DyB 2 C 2 
E| etc.) with quadrupole phase transitions show that 
the magnetic ordering occurs at lower T. A transition 
to an ordered magnetic phase in these rare-earth inter- 
metallic compounds indicates that the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction between localized 
magnetic moments prevails over the the Kondo transi- 
tion mechanism (for a review, see j45|). In zero mag- 
netic field the magnetic susceptibility shows no anomaly 
at Tq for TmTe , CeAg |3^| and a very small anomaly 
for DyB 2 C 2 @, TmAu 2 Within the present work 
we come to the conclusion that quadrupolar order, as 
an electronic charge degrees of freedom driven process 
on one hand, and magnetic properties (Kondo anomaly, 
magnetic order) on the other hand are related indirectly 
via their coupling to lattice displacements. An open ques- 
tion is a possible relation between quadrupolar order and 
hybridization. In our opinion a microscopic derivation 
of the Anderson hybridization Hamiltonian which should 
include details about the symmetry of the lattice site and 
the electronic orbitals constitutes a challenge for further 
work. 
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VII. ERRATUM FOR REFERENCE [21] 



APPENDIX A: 

The Fourier transforms of the intersite electronic 
quadrupole-quadrupole interactions are 3x3 matrices 
where the rows and columns are labeled by the indices of 
the functions S^, k = 1, 2, 3. The structure of these 
matrices depends on the symmetry of the lattice and of 
the Tig functions. The magnitude of the elements de- 
pends on the nature of the electrons (localized 4/, or con- 
duction 6s, 5d) and we abbreviate the electronic indices 
/ or s, d by the label A, B for F = (//) or L = (hh) 



writing v£%, for vffy etc.. We consider elements 



,,ab 

] k k 



■4m 



,iq-X(h) 



(Al) 



Here v AB (h) refers to the elements in real space with 
h = n' — n. Taking the site ft as origin on a fee lattice, 
the index n' (or h) la be ls the twelv e neig hbors. Compare 
with expressions J2.6| ), ( 3.10a ) and ( 3.17 ) for the case A 



{Tig,k), A' = (T 2g ,/c'). Performing the lattice sums and 
using the symmetries of the elements v AB (h) we obtain 



(A2) 



-P AB s^ -p A 



(C X y -\-Cy Z ) 
'Sy. ~f AB C a 



-f3 AB S zx 

-t3 AB Syz 



where CV, = cos(qia/2) cos(qjd/2), and Sij = sm(qia/2) 
x sin(g.,a/2). Here i,j stands for the Cartesian indices 
x, y, z, and a is the cubic lattice constant. Additional 
information on coupling matrices can be found in Ref. 
J47]] where a problem with similar symmetries was con- 
sidered for a molecular crystal. The quantities j AB , 
a AB and (3 AB stand for the matrix elements v AB , (h) for 
h = a(l/2, 1/2,0) with (k = 3, k' = 3), (k = 1, k' = 1) 
and (fc = 1, k 1 = 2), respectively. The interaction ma- 
trix v AB (q) has the largest negative twofold degenerate 
eigenvalue at the Y-point of the Brillouin zone. For in- 
stance for q* = (2?r/a) (1,0,0) X AB = -4 7 AB , where 

V AB (h) — V AB (h) > 0. Hence the quadrupo- 
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AB 



lar interaction parameters X A + are completely specified, 

. 5 

their numerical values are given in Table 2, fourth col- 
umn. 



In Table 3 —/3 (last column, third row) should be re- 
placed with +j3. In the last column of Table 5 —A (first 
row) should be replaced with +A, —fx (second row) - with 
v and — v (third row) - with /i. In Table 2 i — 7, 8 corre- 
spond to Tj, 2 (third row, second column) while i = 13, 14 
- to Tq (fifth row, second column). In equation (4.7d) 
qp(hl) c^\li) has to be replaced by c£\il) c^\li) , on 
the right hand side of equation (5.11) in front of summa- 
tion it should stay i, not |. 



APPENDIX B: 



We study the symmetry properties of the Bloch func- 
tions and rewrite Eq. (|3.]| ) as 

(R\k,a) = ^jR) 

E^'^E^")^-^')). (Bl) 

It 



N *-z 
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Here <fii T (r) — TZi(r)S[(Sl) where S z T (fi) are basis func- 
tions of the irreducible representations of Oh- Here again 
t = (r, fx, k). We recall that 



Y/"(n)=£/?rsf(n), 



(B2a) 



where fii is a unitary matrix. It follows that 

llT (k,a)=Y / Jim (A, a)A mT (B2b) 



We observe that the expression ( 3.10b ) becomes 



Ck'hh(k,a; p,(3) 

= ^2 ^hTi( k : a ) 7i 2 r 2 (P, 0) Cfe' (IlTl, Z 2 T 2 ), (B3c 



where 

cfc' (im, ; 2 r 2 ) = j dn s% (o) s fe , (n) s£ (o). (B3b) 

The Bloch functions can be rewritten as basis functions 
of a d-dimensional irreducible representation of the space 
group Fm3m which is induced by a small representa- 
tion ri p) of the little group of a chosen wave vector 
fei 

k\. Here p labels the small representation. One has 
d = q • t, where q is the number of arms of the star 
of k%: *k% — {ki...ki...k q }, and where t is the dimension 
of the small representation. (If k\ is not a high sym- 
metry vector of the Brillouin zone, q = 48, the order of 
the group Oh, t = 1 and d — 48.) The basis functions 
(R\ki, s p ,a), where s p — 1 — t and fcj € *k\, all have the 
same energy. We can use the small representation index 
p to label the energies. Hence we will write for the band 
index a the composite index (p, v) , where v labels the 
small representations p that occur more than once. We 
consider (see also Appendix B of Ref. I) 

- q - 1 



t q 



s p — 1 i— 1 

for fixed p and fcj S By the generalized Unsold 

theorem p9| this expression is an invariant of Fm3m 
for any r^i. This means that the site symmetry is Oh 
and the density J^k $"(i>k) (r/tfc) (%) sta y s invari- 

ant for any representation T and index n of Oh, (Unsold 
theorem). We therefore obtain the condition 

t q 



5 P7 a ) Ji'T'(ki,s p , ( 



Sp — 1 i — 1 



Sw5 TT i / fcia; r,/j , 



where fci belongs to a basis domain $ [B5| and where / 
is a function of ki and depends on the labels T, /i, a. 

Now we calculat e the single particle expectation values 
{(P l k h (n)) 2 ) , Eq. ( |5.20|) , for conduction elect rons. From 
Eqs. (3.12a,b) we have 



9 33 ft 

x c khi 2 (Pi0 P- q,a)c k >i'j> 2 (h + q,a'-, h,0'). 



(B6) 



Here summation is understood over the repeated band in- 
dices a, a' , (3, P' . The thermal expectation value is calcu- 
lated with the single electron Hamiltonian Uq, Eq. 
We use standard Green's function (GF) techniques |50f| 
and evaluate first the susceptibility (retarded GF) in the 
static frequency limit (z=0) 



X (z = 0) = -{{a}_ ?a a pp - a[ +ifa a Rf3 ,}} z = 



V/3 



E(p-q,a)-E(p,p) 



(B7) 



Here n^ a is the Fermi distribution function. Therefrom 



we get the identity 



V/3 — (1 — n pp) n p-qa 



1 



^E(p-q, a )-E(p,p))/T 



(B8) 



We observe that in the high temperature case the static 
GF {(A; B)) z= q is related to the thermal expectation value 
(AB) by 

(AB) = -T((A:B)) Z=0 . (B9) 
Hence, for \E(p — q, a) — E(p, 0)\ < T, we obtain 

(Pk h (n)Pk' 2 (")>o = EX^ 1 ~ n P^ n P-? a 



a pq 



x c thi 2 (p,P'i p- <f> «) Cfe'/;;^ (p. # P- ( B1 °) 

We next show that symmetry implies that k — k' and 



= l[, h = Expression (B1C) then reduces to (5.20) 
We rewrite expression ( |B10D as 



(Pi lh {n) P P (n))o 4ED 1 " ^) »K« 

x Ckhh&P'-, hi^Ck'i'^&P; h,a). (Bll) 



Using Eq. (B5), we now exploit the fact that the summa- 
tion over the Brillouin zone can be decomposed by dis- 
tinguishing the summation over the basic domain from 
the summation over the components of irreducible repre- 



ss) sentations. Expression (Bll) can be rewritten as 
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(ti h (fi)p l k \ l *(n)) = ((p ( fe l(2 (n)) 2 )o 6 hVl 6 hl , 2 6 kk ,, (B12) [1] 



where 



((p^(n)) 2 ) = EE EEEEE 

a P hiE^PiG* s h 8 s p 3 
i2 

Ckhi 2 (hi,s h ,a; Pj,s p ,f3)\ P{h x a; pi/3). (B13) 



Here we have taken into account the property ( |B5| ) . The 
quantity P is a function of temperature due to the pres- 
ence of the Fermi distribution. We next show that expres- 
sion ( B13 ) is independent of the index of the quadrupolar 
component k (i.e. has the same value for fc=l,2,3.) We 
consider the summation over the components of a pair of 
irreducible representations: 



X = P(h ia - p x p) 

x EEEE |cfcW2(^> Pj,s p ,f3) 

Sh i s p j 



(B14) 



Inserting the explicit expression for Cki 1 i 2 we obtain 

^ = E E P '( h ii a > r i>Mi; Pi, ft T 2 ,M2) 
ri r 2 mi 

xJ2l c t(hn;hT 2 )} 2 , (B15) 



where r, = (ri,/ij,Sj), i = 1,2. We have used relations 
(B3a,b) and (B5). The quantity P' is defined by 



P'(hi,a, Tx,(ii; pi, [3, T 2 ,^ 2 ) 

= P(hia; pi/3) f hia;Ti,Lii /[pi ft r 2 , M2] • (B16) 

Finally we observe that as a consequence of the Wigner- 
Eckart theorem 

^ c* k (l 1 T 1 ;l 2 T2)ck'(liT 1 ;l 2 T2) = Skk'C(r 1 n 1 ;T 2 ^2), 

(B17) 

where C is independent of fc (as a square of a reduced 
matrix element). This proves that X and correspond- 
ingly ((p l k l12 (n)) 2 )o is independent of k. Finally we quote 
a property of 7; r that follows from the completeness of 
the basis functions 

t q 

^2 E E ^2"fiT(kt,s p ,a)~fi> T ,(k l ,Sp,a) = § n ,8 rr ,N. 

(B18) 
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